Correlations Near the Vulcanization Transition: A Renormalizat ion- Group Approach 
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Correlators describing the vulcanization transition are constructed and explored via a renormal- 
ization group approach. This approach is based on a minimal model that accounts for the thermal 
motion of constituents and the quenched random constraints imposed on their motion by crosslinks. 
Critical exponents associated with the correlators are obtained near six dimensions, and found to 
equal those governing analogous entities in percolation theory. Some expectations for how the 
vulcanization transition is realized in two dimensions, developed with H. E. Castillo, are discussed. 
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Introduction. The vulcanization transition (VT) is the 
equilibrium phase transition from a liquid state of matter 
to an amorphous solid state. It occurs when a sufficient 
density of permanent random constraints (e.g. chemical 
crosslinks) — the quenched randomness — are introduced 
to connect the constituents (e.g. macromolecules) , whose 
locations are the thermally fluctuating variables. Whilst 
a rather detailed description of the VT has emerged over 
the past few years within the context of a mean-field ap- 
proximation [|]J|], the picture of this transition beyond 
the mean-field (MF) level is less certain. 

The true (i.e. beyond-MF) critical properties of the 
VT (and the related chemical-gelation transition) have 
mostly been studied through approaches based on gela- 
tion/percolation perspectives [^|-(5), in which the VT is 
directly or indirectly identified with the percolation tran- 
sition. These approaches are associated with one single 
ensemble, accounting either for the quenched disorder or 
for the equilibrium thermal configurations (whose change 
in character mark the transition) . Given that an essential 
aspect of the VT is the impact of the quenched random 
constraints on the thermal motion of the constituents, 
approaches based on a single ensemble cannot directly 
account for the effects of both types of fluctuations and 
are, thus, not entirely satisfactory. 

The purpose of the present Letter is to apply 
renormalization-group (RG) ideas to a model in which 
both the quenched randomness and the thermal fluctua- 
tions of the constituents are naturally and directly incor- 
porated. Our analysis provides a more complete way of 
obtaining the critical exponents, and therefore sheds light 
on the connection between gelation/percolation physics 
and the VT. A detailed account of this work can be found 
in Ref. §. 

Our approach to the VT is based on a minimal Landau- 
Wilson effective Hamiltonian, which was constructed in 
Ref. H and also shown to recover (at the mean-field 
level) the description of both the liquid and emergent 
amorphous solid states, known earlier from the analysis 
of various semi-microscopic models [^,^| . The order pa- 
rameter for the VT is a function that encodes both the 
gel fraction q (i.e. the fraction of monomers localized) 
and the distribution of localization length of the localized 



monomers (as well as other diagnostics). Support for the 
mean-field picture of the solid state has emerged from 
extensive molecular dynamics computer simulations of 
three-dimensional, off-lattice, interacting, macromolecu- 
lar systems, due to Barsky and Plischke JTcj . 
Modeling the VT — the order parameter. The appropri- 
ate (dimensionless) order parameter for the VT, capa- 
ble inter alia of distinguishing between the liquid and 
amorphous solid states of randomly crosslinked macro- 
molecular systems (RCMSs), is the following function of 
n wave- vectors {k 1 , . . . , k™}: 
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where N is the total number of macromolecules, Cj(s) 
(with j = 1, . . . , N and < s < 1) is the position in d- 
dimensional space of the monomer at fractional arclength 
s along the j th macromolecule, (• • -) x denotes a thermal 
average for a particular realization \ of the quenched dis- 
order (i.e. the crosslinking), and [• • •] represents a suit- 
able averaging over this quenched disorder. It is worth 
emphasizing that the disorder resides in the specification 
of what monomers are crosslinked together: the result- 
ing constraints do not explicitly break the translational 
symmetry of the system. 

Modeling the VT — the minimal model. Following Deam 
and Edwards , we adopt a constraint distribution ap- 
propriate to the situation of instantaneous crosslinking of 
the equilibrium melt or solution, and invoke the replica 
trick to incorporate the consequences of the permanent 
random constraints. Thus, we are led to the need to work 
with the n — > limit of systems of n+ 1 replicas. The ad- 
ditional (i.e. zeroth) replica incorporates the constraint 
distribution. 

The form of the minimal model can be determined, in 
the spirit of the Landau approach, from the nature of 
the order parameter and the symmetries of the effective 
(i.e. pure but replicated) Hamiltonian, along with the as- 
sumptions of the analyticity of this Hamiltonian and the 
continuity of the transition Q . This scheme leads to the 
following minimal model, which takes the form of a cu- 
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bic field theory involving the order parameter field £l(k) 
living on (n + 1 )-fold replicated d-dimensional space: 

/ cx lim n- 1 In \Z n \ , [Z n ] oc / P t fiexp(-5„) , (2a) 

TWO J 

s n ({n}) = N J2 (-«T + ^|fc| 2 )|fi(fc)| 2 

A; 6 HRS 

-Ng Yl n (*i) n (*a)n(fc3)^ l+fe+M . (2b) 

Here, r is the control parameter, which measures the 
reduced density of random constraints, and the coeffi- 
cients a, b and g encode the microscopic details of the 
system. We denote averages weighted with exp(— <S„) by 
(■ • •} <s , use the symbol k to denote the replicated wave- 
vector {k°, k 1 , . . . , k™}, and define the extended scalar 
product k ■ c by k° ■ c° + k 1 • c 1 + • • ■ + k™ • c™. The sym- 
bol k £ HRS indicates that a summation over replicated 
wave- vectors is restricted to those containing at least two 
nonzero component- vectors k". (We say that this kind 
of wave- vector lies in higher-replica-sector, i.e. the HRS.) 
This condition reflects the important fact that no crys- 
talline order (or any other kind of macroscopic inhomo- 
geneity) is present in the vicinity of the VT; it changes 
the symmetry of the effective Hamiltonian and, hence, 
can (and in fact does) play a crucial role. Consequently, 
the functional integral is to be performed over fields obey- 
ing the linear constraint of lying in the HRS. 
Ginzburg criterion. To estimate the range of r about the 
critical value, within which the effects of order-parameter 
fluctuations are relatively strong, we follow the conven- 
tional strategy of constructing a loop expansion (in the 
present setting, an expansion in the inverse monomer 
density) for the 2-point vertex function to one-loop or- 
der, examining its low- wave- vector limit, thus obtaining 
the inverse susceptibility S. 




FIG. I. One-loop correction to the 2-point vertex function. 

The only one-loop correction comes from the diagram 
shown in Fig. []]; thus, we find 

^-- 2 ^"^(^ 18 ^ ( " 25T)( " 6)/2 )' 

where r c > is the critical value of r, shifted due to 
the inclusion of one-loop corrections, and Jd is an unim- 
portant dimensionless number. Equation (^|) shows that 
the upper critical dimension for the VT is six. To deter- 
mine the physical content of the Ginzburg criterion, we 
invoke the values of the coefficients appropriate for the 



semi-microscopic model of RCMSs, and find the follow- 
ing form of the Ginzburg criterion: for d < 6, fluctuations 
cannot be neglected for values of r satisfying: 

\T\<{L/ir^^/g*y 2/ ^ d \ (4) 

Here, L is the arclength of each macromolecule, I the per- 
sistence length, and ip = (N/V)(L/£)£ d the volume frac- 
tion. Equation (|J) shows that the fluctuation-dominated 
regime is narrower for longer macromolecules and higher 
densities (for 2 < d < 6). Such dependence on the de- 
gree of polymerization L/£ is precisely that argued for by 
de Gennes on the basis of a percolation picture jl2| . 
Epsilon expansion. We now apply the e (= 6 — d) expan- 
sion to derive the RG flow equations near the upper criti- 
cal dimension, and discuss the resulting fixed-point struc- 
ture and critical exponents. In order to streamline the 
presentation, we may, by suitably redefining the scales of 
Q(k) and k in Eq. (pb|), absorb the coefficients a and b 
(i.e. set a = 1, b = 1). 




FIG. 2. Contributing one-loop diagrams. Full lines 
indicate bare HRS correlators for short-wavelength fields 
(i.e. fields lying in the momentum shell); wavy lines indicate 
long-wavelength fields. 

We shall be working to one-loop order and, correspond- 
ingly, the contributing diagrams are those depicted in 
Figs. ||a,b. After taking proper care of the constraints on 
the wave-vector summations, and passing to the n — > 
limit, the flow equations are (with higher order terms 
omitted) found to be 

dr/dlnb = 2r - C g 2 - C' q t g 2 - C x t g 1 , (5a) 
dg/d\nb = g(e/2-C 3 g 2 -^C ig 2 ), (5b) 

dz/d\nb= ^(d + 2-dg 2 ), (5c) 

where b is the length-rescaling factor, z is the field- 
rescaling factor, and the (constant) coefficients in 
the flow equations are given by (Co, Cq, C\, C3) = 
(V r /iV)(S' 6 /(27r) 6 )(9A 2 , 36, -6, 72) in which S e is the sur- 
face area of a 6-dimensional sphere of unit radius and A 
is the cut-off for replicated wave-vectors. 

Proceeding in the standard way, we find that: (i) For 
e negative (i.e. d > 6), there is only the Gaussian 
fixed point (GFP), (t*,^) = (0,0), at which the 
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exponents take on their classical value: v = 2, 
rj = 0. (ii) For e positive (i.e. d < 6), in ad- 
dition to GFP, a new fixed point — the Wilson-Fisher 
fixed point (WFFP) — emerges, located at (r*,g 2 ) = 
((A 2 /28),(l/126)((27r)7S* 6 )(V r /Af)- 1 )£, controlling the 
critical behavior. The resulting critical exponents are 
(to first order in e) v^ 1 = 2 — (5e/21) and rj = — e/21. 

A standard scaling argument leads to the value of the 
critical exponent j3 for the gel fraction: for d > 6 we find 
(3 = 1; for d < 6 we find /3 = u(d + 2 - rj)/2 = 1 - e/7. 
In fact, under the (not unreasonable) assumption that, 
in the ordered state, the fluctuation correlation length- 
scale is the same as the localization lengthscale, we can 
go further and propose a general scaling hypothesis for 
the (singular part of) (n(k)) s , viz., 
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which is supported by the mean-field result M. There 
are, however, some fascinating subtleties here, which may 
instead yield multi-fractal characteristics for vulcanized 
matter similar to those explored by Harris and Luben- 
sky in the setting of random xy models and resistor net- 
works |]|. 

Order-parameter correlator and its physical significance. 
Both above and below six dimensions, the critical expo- 
nents u, rj and (3 of the VT, computed above, have iden- 
tical values (at least to first order in e) to their physical 
counterparts in percolation theory (as computed, e.g., via 
the Potts field theory |l4|]). To establish this correspon- 
dence, we now explore the physical content of the order 
parameter correlator (We remark that the identifi- 
cation of the VT (3 with the j3 of percolation is evident.) 
To proceed, let us consider the construct 
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where the term in the second line that is delimited by 
square brackets is the microscopic counterpart of the cor- 
relator of a typical component of the order parameter, 
i.e., (n(0, . . . , 0, k, t)* 0(0, . . . , 0, k, t)) S . The first ther- 
mal average in Eq. ((7a|) accounts for the likelihood that 
monomers (j, s) and (j',s') are respectively to be found 
around r and r'; the second describes the correlation be- 
tween the respective fluctuations of monomer (j, s) about 
r and monomer (J', s') about r'. 

The small-t limit of Ct(r — r') addresses the con- 
nectedness of clusters of mutually crosslinked macro- 
molecules 1 16 . To substantiate this claim, we examine 
Eq. (|7a|), and consider the contribution to Ct(r — r') from 
pairs of monomers that are in the same cluster and pairs 
that are in different clusters. (We assume that the sys- 
tem has only short-range interactions.) For a generic pair 
of monomers that are in the same cluster, we expect that 



(expit- (cj(s) — Cji(s'))) 



1, and that (for |r-r'| < R s ) 
j («'))) ~ V- 1 R- d . Then the 
total miro-cluster contribution to Ct(r — r') is at least 
of order (N/V) 2 R~ d . On the other hand, a similar anal- 
ysis shows that the total m£er-cluster contribution is at 
most of order (N/V^V -1 . Therefore, the intra-cluster 
contribution dominates Ct(r — r') in the thermodynamic 
limit. (This is also true in the case of generic t, and thus 
we can always ignore the inter-cluster contribution.) In 
other words, in the small-t limit, a pair of monomers lo- 
cated at r and r' contribute unity to Ct(r — r') if they 
are on the same cluster and contribute zero otherwise. 



This limit of Ct(r — r') plays the same role as the pair- 
connectedness function defined in (the on-lattice version 
of) percolation theory [||. [For t = 0, Ct(r — r') is sim- 
ply {V/N times) the real-space density-density correla- 
tor, and is not of central relevance at the VT.] 

In the case of general t, Ct(r — r') addresses the ques- 
tion: If a monomer near r is localized on the scale i _1 
(or more strongly), how likely is a monomer near r' to 
be localized on the same scale (or more strongly)? From 
Eq. ([[a]) we see that t~ Y serves as a cutoff to the range of 
the correlator, so that all pairs of monomers that are rel- 
atively localized on a length-scale much larger than 1/t 
do not contribute to Ct (r — r' ) . 

Given Ct(r — r'), we can build a divergent susceptibil- 
ity by integrating Ct(r — r') over space and passing to 
the t — > limit: 



f d d r d d r' , 
N lim (n(0, ..,0,t, -t) 

n— >0 



(8) 



0(0,..,0,t,-t)) 6 ~(-r)" 



This quantity is measure of the spatial extent over which 
pairs of monomers are relatively localized, no matter how 
weakly, and thus diverges at the VT. 
Concluding remarks. As we have seen, the VT expo- 
nents obtained by the e expansion turn out to be numer- 
ically equal to those characterizing physically analogous 
quantities in percolation theory, at least to first order in 
e. This equality between exponents seems reasonable, 
in view of the intimate relationship between percolation 
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theory and the connectivity of the system of crosslinked 
macromolecules. Such a connection has long been recog- 
nized, and supports the use of percolative approaches as 
models of certain aspects of the VT . 

A convenient formulation of the percolative approach 
is via the Potts model in its one-state limit [|17|. It is 
therefore worth considering similarities and differences 
between the minimal model of VT, Eq. (2b), and the 
minimal field theory for the Potts model. Both have a 
cubic interaction and both involve an in limit. De- 
spite the similarities, however, there exist many differ- 
ences: (i) The Potts field theory has a multiplet of n real 
fields on d-dimensional space; the VT field theory has 
a real singlet field living on (n + l)c?-dimensional space, 
(ii) The Potts field theory emerges from a setting in- 
volving a single ensemble, whereas the VT field theory 
describes a physical problem in which two distinct en- 
sembles (thermal and disorder) play essential roles. As 
such, the vulcanization field theory is a more natural and 
direct approach, in the sense that it is capable of provid- 
ing a unified theory of both the transition and the struc- 
ture of the solid state, (iii) The symmetry structures 
are different; most strikingly, the nature of the spon- 
taneous symmetry breaking is different. The percola- 
tion transition (in its Potts representation) involves the 
spontaneous breaking of the (n — > limit) of a discrete 
(n + I)-fold permutation symmetry. By contrast, the VT 
involves the spontaneous breaking of the (n — > limit 
of the) continuous symmetry of relative translations and 
rotations of the n + 1 replicas; the permutation symme- 
try remains intact in the amorphous solid state, as does 
the symmetry of common translations and rotations of 
replicated space. The fact that the critical exponents 
near the upper critical dimension are the same for both 
theories (at least to order e) despite these apparent dis- 
tinctions can be reconciled by delineating between three 
logically distinct physical properties pertaining to amor- 
phous solidification transition: (i) macroscopic network 
formation; (ii) changes in the nature of thermal motion — 
i.e. random localization; and (iii) the acquisition of rigid- 
ity. Within mean-field theory (and hence above six spa- 
tial dimensions), these three properties go hand- in-hand, 
emerging simultaneously at the VT. At and below six 
dimensions they appear to continue to go hand-in-hand 
(although we have not yet investigated the issue of rigid- 
ity beyond mean-field theory) until one reaches two di- 
mensions, where (as we shall discuss shortly) we believe 
this broad picture changes. Thus, in higher dimensions it 
appears that when both ensembles are incorporated the 
disorder fluctuations play a more important role than the 
thermal fluctuations. 

This brings up the interesting issue of the nature of the 
VT at the neighborhood of two dimensions — the lower 
critical dimension of the VT [lq| . (The ideas reported in 
this paragraph result from an ongoing collaboration with 
H. E. Castillo Indeed, there is a conventional perco- 

lation transition in two dimensions, whereas the thermal 
fluctuations are expected to destroy the random localiza- 



tion. It is tempting to speculate |19| that in two dimen- 
sions an anomalous type of VT happens simultaneously 
with percolation transition: As the constraint-density is 
tuned from below to above its critical value, the order 
parameter would remain zero, whereas its correlations 
would change from decaying exponentially to decaying 
algebraically: one would have a quasi-amorphous solid 
state — the random analog of a two-dimensional solid |20| . 
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